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Abstract 

We report detailed measurements of the interlayer magnetoresistance of the 
layered organic superconductor k-(BEDT-TTF)2Cu(SCN)2 for temperatures 
down to 0.5 K and fields up to 30 tesla. The upper critical field is deter- 
mined from the resistive transition for a wide range of temperatures and field 
directions. For magnetic fields parallel to the layers, the upper critical field 
increases approximately linearly with decreasing temperature. The upper crit- 
ical field at low temperatures is compared to the Pauli paramagnetic limit, at 
which singlet superconductivity should be destroyed by the Zeeman splitting 
of the electron spins. The measured value is comparable to a value for the 
paramagnetic limit calculated from thermodynamic quantities but exceeds 
the limit calculated from BCS theory. The angular dependence of the upper 
critical field shows a cusp-like feature for fields close to the layers, consistent 
with decoupled layers. 
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I. INTRODUCTION 



The layered organic molecular crystals k-(BEDT-TTF)2X (where BEDT-TTF is 
bis-(ethylenedithia-tetrathiafulvalene) and X is an anion (e.g., X=l3, Cu[N(CN)2]Br, 
Cu(SCN)2)) are particularly interesting because they are strongly correlated electron systems 
with similarities to the high-Tc cuprate superconductors including unconventional metallic 
properties and competition between antiferromagnetism and superconductivity |[T|-Q. Fur- 
thermore, they are available in high purity single crystals and, in contrast to the cuprates, 
their lower superconducting transition temperature (Tc ~ lOK) makes experimentally acces- 
sible in steady magnetic fields properties such as the upper critical field and Shubnikov-de 
Haas oscillations 

Recently it has been argued that a minimal theoretical model that can describe these 
materials is a Hubbard model on an anisotropic triangular lattice with one hole per site 
PJ^. Calculations at the level of the random-phase approximation [0] and the fiuctuation- 
exchange approximation suggest that at the boundary of the antiferromagnetic phase this 
model exhibits superconductivity mediated by spin fiuctuations. As the anisotropy of the 
intersite hopping varies the model changes from the square lattice to the isotropic triangular 
lattice to decoupled chains 0. The wavevector associated with the antiferromagnetic spin 
fiuctuations changes [^] and the superconductivity has been predicted to change from d-wave 
singlet (as in the cuprates) to s-wave triplet in the odd-frequency channel 0. 

Experimental results that are consistent with unconventional superconductivity include 
the temperature dependence of the nmr relaxation rate 1/Ti (including the absence of a 
Hebel-Slichter peak) |lTO| , p!T| , the temperature and magnetic field dependence of the electronic 
specific heat 0], the temperature dependence of the thermal conductivity, and the 
sensitivity of Tc to disorder . 

The temperature dependence of the nmr Knight shift (which measures the electron spin 
susceptibility) in the superconducting state provides a means to distinguish triplet and 
singlet pairing. For triplet pairing the Knight shift does not change on entering the super- 
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conducting state, whereas for singlet pairing the Knight shift goes to zero as the temperature 
decreases to zero. The Knight shift of ^'^C nmr on the X=Cu[N(CN)2]Br is consistent with 
the latter. In contrast, the Knight shift of ^''O nmr on Sr2Ru04 is consistent with the former. 
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If the superconductivity is spin singlet then the upper critical field cannot exceed the 
paramagnetic limit Hp, also known as the Pauli limit or Clogston-Shandrasekhar limit 
T6| , [T7[| . Above Hp the Cooper pairs are destroyed by the Zeeman splitting produced by the 



magnetic field coupling to the electronic spins. For weak-coupling BCS theory 

Hp = H^^' r:^^^^. (1) 

For Tc = 10 K, as in the material studied here, this gives Hp'~^^ = 18 T. Strong coupling 
effects [0 and d-wave pairing |jl9[ only change this value of Hp slightly. In most supercon- 



ductors the paramagnetic limit is irrelevant because the superconductivity is destroyed at 
much lower fields due to the frustration of the orbital degrees of freedom associated with the 
formation of vortices. However, in layered superconductors with fields parallel to the layers 
the vortices can fit between the layers and paramagnetic limiting can become important. 

m 

Previous determinations of the upper critical field of the fi;-(BEDT-TTF)2X family 
|2l|-p6| have mostly focussed on measurements of the slope "^^^t^^ near T^. The values 
obtained for X=Cu[N(CN)2]Br and X=Cu(SCN)2 are in the range 10 to 20 T/K. Using 
the WHH formula for a three-dimensional superconductor, this very large slope would 



suggest a zero-temperature Hc2{T = 0) = 0.7Tc ^^2t^^ = 70-140 T, which is well above the 
BCS Pauli limit. A previous transport measurement on the X=Cu(SCN)2 salt was carried 
out in pulsed magnet fields. A quasi-linear temperature dependence was found with 
Hc2 ~25 T and the authors concluded that the upper critical field exceeded the Pauli limit. 
A study of the upper critical field of X=Cu(CN)[N(CN)2] determined from the resistive 
transition found an upper critical field of about 25 T for fields parallel to the layers. Studies 
on the lower Tc organic compounds such as the K-(BEDT-TTF)2l3 |^0[|, /3-(BEDT-TTF)2l3, 



and /3-(BEDT-TTF)2lBr2 [^] have found that the Hc2 at zero temperature hes below or 



close to the Pauli paramagnetic limit predicted by BCS theory. Similar paramagnetic field 
limited ifc2 have been reported in the cuprate YBa2Cu307_5 and the heavy fermion 
superconductors UPd2Al3 [p3| . 

If there is paramagnetic limiting there is theoretically the possibility that as the magnetic 
field is increased at low temperatures there is a first-order phase transition into non-uniform 
superconducting state, originally proposed by Fulde and Ferrell and Larkin and Ovchinikov 
|3^ . As the dimensionality of the system decreases the magnetic field range over which 



this phase is stable increases. [35| Such a first-order phase transition was recently seen in 



ultrathin beryllium films. |36] It is still controversial about whether this phase does exist 



in UPd2Al3 [Q. On the other hand, if the superconductivity is triplet there is also the 



possibility of re-entrant superconductivity at high fields such that Tf.{H) actually increases 
with increasing field p7| , p8| . 

In this paper we report the measurement of the interlayer resistivity of X=Cu(SCN)2 
down to 0.5 K and up to 30 T for a range of field directions. For magnetic fields parallel to the 
layers, the upper critical field increases approximately linearly with decreasing temperature 
to values that clearly exceed the BCS Pauli limiting field (|Tp, but are consistent with the 
paramagnetic limit, estimated directly from the superconducting condensation energy. The 
upper critical field as a function of angle shows a sharp cusp for fields almost parallel to 
the layers, consistent with two-dimensional decoupled layers. We find no evidence of a first 
order phase transition as a function of field at low temperatures. 

II. THEORETICAL BACKGROUND 

We now briefiy summaries some theoretical results concerning the upper critical field 
which we will use later in interpreting our results. A more complete discussion can be found 



in Reference p9| 
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A. Angular dependence of the upper critical field 



Anisotropic Ginzburg-Landau theory is valid when the coherence length perpendicular 
to the layers, is much larger than the interlayer spacing. It predicts that the dependence 
of the upper critical field on the angle 6 between the field and the normal to the layers is 



H. 



c2± 



+ 



H, 



c2\\ 



1, 



(2) 



where Hc2± and Hc2\\ are the upper critical field for fields perpendicular and parallel to the 
layers, respectively. The perpendicular upper critical field is determined by ^|[, the coherence 
length parallel to the layers, 

$0 



c2± 



(3) 



where $o is the fiux quantum. The coherence lengths parallel and perpendicular to the 
layers are related by 

Hc2± 



H, 



(4) 
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Klemm, Luther, and Beasley considered the upper critical field of layered supercon- 



ductors when the layers were infinitely thin. |20| For both Lawrence-Doniach theory and 
microscopic theory, they found that for fields parallel to the layers, if the interlayer coupling 
is sufficiently weak the upper critical field diverges at low temperatures unless spin-orbit 
effects or paramagnetic limiting is present. This is because the Josephson vortices associ- 
ated with the field parallel to the layers have no normal core and can fit between the layers. 
Bulaevskii and Schneider and Schmidt considered a more general model where the 
layers have a finite thickness d, resulting in a finite upper critical field 

Hc2± d 



H, 



c2\\ 



12^11 



(5) 



They also found that if the coupling between the layers is sufficiently weak then the angular 
dependence of the upper critical field is given by 



H,2{9) cos(g) 

Hc2± 

This same angular dependence was found earlier for thin two-dimensional films by Tinkham 
using a simple fluxoid quantization argument. The main difference from the anisotropic 
three-dimensional result is that at 9 =90°, Hc2{9) from Eqn. (Q) is smooth or bell-shaped 
with '^^fg^^ = whereas Hc2{9) from Eqn. (^) has a cusp at 9 =90°. 

If the upper critical field is determined solely by coupling of the field to the spins then 
it will be independent of the field direction. Bulaevskii considered the case where the 
paramagnetic limit is larger than the upper critical field for fields perpendicular to the layers 
but smaller than the upper critical field determined by orbital effects for fields parallel to 
the layers. The angular dependence is then given by 



Hc2{9) cos{9) 
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where Hc2\\ = Hp. This also results in an Hc2 versus 9 curve which has a cusp at 9 =90°. 
Indeed the angular dependence is difficult to distinguish from Eqn. (H). 
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B. Estimating the paramagnetic limiting field 

The metallic phase has a finite Pauli spin susceptibility Xe compared to the vanishing 
susceptibility (at zero temperature) of a spin singlet superconducting state. Hence, it will 
be energetically favorable to destroy the superconducting state when the magnetic energy 
density gained by the difference in susceptibilities exceeds the superconducting condensation 



energy density Uc- The critical field Hp at which this occurs is given by [16 



U. = ^XeHl (8) 

where /Xq is the magnetic permeability of free space. 

In BCS theory the condensation energy density is Uc = ^N{Ef)A{OY where N{Ef) 
is the metallic density of states and A(0) = l.TQkpTc is the zero-temperature energy gap. 
Making using of these relations and Xe = {f^B)'^N{EF) we obtain the expression (|l|) for Hp. 
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Many-body effects. In the k-(BEDT-TTF)2X crystals there are significant many-body 
effects; the electron effective mass m* determined from magnetic oscillations can be two to 
five times larger than that predicted by band structure calculations. [EIB| The effect of this 



on the paramagnetic limit needs to be taken into account. Perez- Gonzalez finds that 
the paramagnetic limiting field is enhanced by a factor of m*/mh. However, he did not take 
into account the simultaneous effect on the Zeeman splitting: the g factor changes to g*. 
When this is done one finds that within a Fermi liquid framework the Pauli limit is actually 
reduced from (0) by a factor of g*/g. |^ This ratio can be estimated from thermodynamic 
measurements or from the spin-splitting of magnetic oscillations. [0] The values obtained 
by these two methods for X=Cu(SCN)2, are 0.8 and 1.4, respectively. 



Alternatively, we can make a theory-independent estimate of Hp by using (Q) and the 
experimentally determined condensation energy density and spin susceptibility. This method 
of determining Hp is very attractive because it does include all the many-body effects 
(without assuming a Fermi liquid picture) and does not assume the validity of any particular 



theory of superconductivity for the material in question. Haddon et al. |^ found Xe = 
4.3 X 10~^ emu per mole (corresponding to a density of states of 7 states per (eV molecule)) for 
the X=Cu(SCN)2 salt. By a reanalysis of Graebner et al.'s [^] specific heat data Wosnitza 
1^ evaluated the condensation energy density in terms of the thermodynamic critical field 
Bth = 90 mT where Uc = '^B'^h- Taking the unit cell volume of 1695 and two (BEDT- 
TTF)2X units in each unit cell gives Bp = 30 ± 5 T. The uncertainty is estimated based on 
the uncertainty in the values for the condensation energy and the susceptibility. 

III. EXPERIMENTAL DETAILS 

Single crystals of /t-(BEDT-TTF)2Cu(SCN)2 were synthesized by the electrocrystalliza- 
tion technique described elsewhere |T^. The interlayer resistance was measured with use of 



the four probe technique. Contact of the gold wires to the sample was made with a Dupont 
conducting paste or graphite paste. Typical contact resistances between the gold wire and 
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the sample were about 10 Q,. A current of 1 iJ,A was used to ensure linear /- F characteristics. 
The voltage was detected with a lock-in amplifier at low frequencies of about 312 Hz. To 
avoid pressure effects due to solidification of grease, the sample was mechanically held by 
thin gold wires. The data presented in this work were taken in a ^He system with field up 
to 30 T at the National High Magnetic Field Laboratory at Tallahassee. The sample can be 
rotated in the field and the orientation was determined by using a Hall probe at low fields. 

IV. RESULTS 

Shown in Fig. 1 is a typical field dependence of the interlaycr resistance plotted in a 
semi-log scale at a temperature of 4.2 K. The field is applied parallel to the planes. The 
resistive transition in parallel field is typical of the low dimensional organic superconductors 
with a broad transition width in field and a large positive magnetoresistance in the normal 
state. The superconducting transition or the upper critical field Hc2 is defined at the 1 fl 
level. To check the validity of this criteria, the critical field will be compared with that 
obtained by a more conventional definition. Shown in the inset are the same data in a 
linear scale. The two lines are extrapolations of the normal state magnetoresistance and the 
superconducting transition with the upper critical field Hc2 * defined at the crossing point of 
the two lines. 

Fig. 2 is an overlay of resistive transitions in parallel field at different temperatures from 
T = 0.5 K to 10.2 K. With increasing temperature, the curves shift to the left toward lower 
critical fields. The transition curves are nearly parallel for all temperatures in the semilog 
scale. Hc2 is almost the mid-transition point as in a conventional superconductor, where 
parallel transitions are seen but in a linear scale. 

The temperature dependences of the two fields Hc2 and Hc2 * are shown in Fig. 3. Within 
the scatter of the points, the two upper critical fields have nearly the same linear temper- 
ature dependence with dHc2/dT 3 TK~^. The offset in the superconducting transition 
temperature is due to the different definitions. The upper critical fields at zero temperature 
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are about 30 T and 33 T for Hc2 and Hc2 *, respectively. The dashed hne is the Pauh hmit 
Hp = 18.4 T, calculated from Eqn. (|I]) with = 10 K. Clearly, Hp defined this way is 
well below the measured upper critical fields at low temperatures. On the other hand, Hc2 
is consistent with our estimate of Hp from thermodynamic quantities. 

To look at the anisotropy of the upper critical field, a systematic measurements have been 
taken as a function of angle 6, defined between the field direction and the normal of the 
plane. Plotted in Fig. 4 is an overlay of resistive transitions as a function of field at different 
angles. The six curves are representative of the angular dependence from field parallel to the 
layers {6 = 90") to normal to the layers {6 = 180°). With increasing 6, the field dependence 
of the resistive transition is drastically changed. At ^ = 91.50°, Hc2 is decreased by about 4 
T. At = 96.64°, a shoulder-hke feature is developed in R(H) with a corresponding decrease 
in Hc2 by about 12 T. The shoulder-like structure develops into a well defined peak at 6 
= 178° with the occurrence of the Shubnikov de-Hass (SdH) oscillation in the resistance 
at high fields. It should be noted that unlike for fields parallel to the layers, the resistive 
transition is relatively insensitive to the angles near 6 = 180°. 

The inset in Fig. 4 shows an expanded view of the resistive transitions at angles close to 
6 =90° direction. With a slight increment in 6, the transition is drastically broadened. The 
field component parallel to planes is almost constant for all angles shown in the inset and 
the maximum out of plane field component is about 0.5 T at = 91.50° and H = 30 T. 

Hc2 defined at the 1 Q level as a function of angle is shown in Fig. 5 at a temperature 
of 1.56 K. Clearly, Hc2 decreases rapidly away from the parallel to the plane direction and 
is nearly saturated above 140°. The three lines are fits to the three-dimensional anisotropic 
model (Eqn. (||)) and the decoupled layer results (Eqn. (||) and (0)). The fits to Eqn. (H) 
and d^) are indistinguishable in the scale shown. While all three fits seem reasonable at 
first sight, clear deviations are seen very close to 6 =90°, as shown in the inset. A cusp-like 
feature is observed experimentally, as in the fit to the decoupled layer model, while the 3D 
fit is rounded with a negative curvature at the top. A better agreement with the data for 
the decoupled layer model at large angles is also evident with the 3D fit lying systematically 
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under the data. At 1.56 K, the 2D fit gives ifc2± = 2.27 T and Hc2\\ =24.5 T. Wanka et al 
also found that the angular dependence for the X=l3 salt was fit best by Eqn. (||). 

V. DISCUSSION 

Our value of Hc2± = 2.3 T at 1.56 K can be compared with the value of about 1.8 T found 
from the irreversibility line deduced from torque measurements. From the perpendicular 
upper critical field value of 2.3 T and Eqn. @ we deduce an intralayer coherence length 
of 120 A. The anisotropic three-dimensional theory (Eqn. (^)) and the measured ratio of 
the upper critical fields gives a perpendicular coherence length of ^± ~ 13A. Since this is 
comparable to the interlayer spacing of 15A we cannot expect the theory to apply. Hence, 
it is not surprising that the angular dependence is not described by Eqn. (||). 

If instead we consider the model of weakly coupled layers and use Eqn. (^) for the ratio 
of the critical fields we deduce that the thickness of the superconducting layer is ci = 40A. 
Clearly, this is unrealistic because it should be smaller than the interlayer spacing. A more 
realistic value would be a few A. This suggests that the parallel upper critical field being 
determined by paramagnetic limiting rather than orbital effects is more realistic. 

Because of the extremely sensitive angular dependence of the resistive transition, a 
shoulder-like feature is developed in the resistive transition a few degrees away from the 
parallel to the plane direction. The upper critical field Hc2* can only be defined close to 
the planes. While the magnitude of Hc2* is larger than Hc2, as expected, it is difficult to 
distinguish the 2D and the 3D models with the available data. Hc2* decreases quasi-linearly 
with angle within the errors. 

The upper critical field determined from transport measurements has been under a lot 
of debate in the cuprate superconductors [^. For field perpendicular to the planes, Hc2(T) 
defined at certain fractional normal state resistance typically gives rise to a positive curva- 
ture at low temperatures. Various mechanisms have been proposed for the unconventional 
temperature dependence. However, it has been suggested that the Hc2 thus defined corre- 
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spends to the irreversibility or vortex melting line. For fields parallel to the layers, a vortex 
moving along the plane encounters negligible pinning as there is no normal core associated 
with Josephson vortices. Magnetization is practically always reversible in this orientation. 
The resistive onset field is clearly well separated from irreversibility field and refiects the 
true upper critical field. 

In the case of Sr2Ru04 and the quasi-one-dimensional organic superconductor 
(TMTSF)2X where X=C104 and PFe [Q, the upper critical field in the plane has been 
found to exceed the Pauli limit, calculated from BCS theory. Combined with the strong de- 
pendence of the transition temperature on the impurity concentration and the temperature 
dependence of the Knight shift, triplet pairing or p-wave has been suggested in these systems. 
However, the quasi-linear temperature dependence observed here for both Hc2 and Hc2* is 
remarkably different from that of Sr2Ru04 and Bechgaard salts. For both Sr2Ru04 and 
(TMTSF)2C104, the H^2 is found to saturate for T/T, < 0.2-0.4. While for (TMTSF)2PF6, 
ifc2(T) along both a and V axes where X=C104 displays a diverging temperature dependence 
near T = K. 



VI. CONCLUSIONS 

In summary, for fields parallel to the layers we have observed an upper critical field de- 
termined from resistive transition which is comparable to the paramagnetic limit estimated 
from thermodynamic quantities but is considerably larger than that calculated from BCS 
theory. There is no evidence of a first order transition in the field dependence of the resis- 
tivity, which would occur if there was a transition to a Fulde-Ferrell phase. The observed 
anisotropy of the upper critical field is much less than would be predicted by a model without 
paramagnetic limiting. The upper critical field determined is quasi-linear with temperature. 
The angular dependence of the resistive transition is consistent with the highly anisotropic 
nature of the title compound with a cusp-like angular dependence for field near the plane. 
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Figure captions 

Fig.l Determination of the upper critical field. The main figure shows the interlayer 
resistance as a function of magnetic field on a semi-logarithmic scale, the upper critical field 
being defined as the field at which the resistance is 1 ohm. In the inset the upper critical 
field Hc2* is determined by linear extrapolation. The temperature is 4.2 K and the field is 
parallel to the layers. 

Fig. 2 The field dependence of the interlayer resistance is shown for various temperatures. 
The field is parallel to the layers. 

Fig. 3 The temperature dependence of the upper critical field for fields parallel to the 
layers. The dashed line marks the Pauli paramagnetic limiting field predicted by BCS 
theory. The two values of the upper critical field correspond to the two different methods of 
determination (see Fig. 1). The solid lines are guides for eyes. 

Fig. 4 Dependence of the resistive transition on the field direction. The field dependence 
of the interlayer resistance is shown at various field directions. The angles given denote the 
angle between the field and the normal to the layers. The temperature is 1.56 K. The inset 
shows how the resistive transition becomes significantly broader as the field is moved slightly 
away from the plane of the layers, for which 9 = 90°. 

Fig. 5 Dependence of the upper critical field on the field direction. The dashed line is 
a fit to the anisotropic three-dimensional Ginzburg-Landau result (Eqn. (^), and the solid 
line to the results for weakly coupled layers. (The curves corresponding to Eqn. (^, (^ 
which neglect and include paramagnetic limiting, respectively are indistinguishable.) The 
inset is an expanded view near 6 = 90°, showing that the weakly coupled layer models give 
the best fit. The temperature is 1.56 K. 
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